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Abstract 

We give further support to Smirnov's conjecture on the exact kink S-matrix for the 
massive Quantum Field Theory describing the integrable perturbation of the c = 0.7 mini- 
mal Conformal Field theory (known to describe the tri-critical Ising model) by the operator 
(j)2,i- This operator has conformal dimensions ^) and is identified with the subleading 
magnetic operator of the tri-critical Ising model. In this paper we apply the Thermody- 
namic Bethe Ansatz (TBA) approach to the kink scattering theory by explicitly utilising 
its relationship with the solvable lattice hard hexagon model. Analytically examining the 
ultraviolet scaling limit we recover the expected central charge c = 0.7 of the tri-critical 
Ising model. We also compare numerical values for the ground state energy of the finite 
size system obtained from the TBA equations with the results obtained by the Truncated 
Conformal Space Approach and Conformal Perturbation Theory. 
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1 Introduction 



Since the work [jT| of A. B. Zamolodchikov, it is known that certain perturbations of conformal 
field theories (CFT's) lead to completely integrable models of massive quantum field theory 
(QFT). In this paper we consider the relevant perturbation of the c = 0.7 minimal CFT 
(known to describe the tri-critical Ising model (TIM) [^) by the operator 02, i- This operator has 
conformal dimensions (A, A) = (^, ^) and is identified with the subleading magnetic operator 
of the TIM. The action of the perturbed tri-critical Ising model (PTIM) can be written as 

ApTIM = AtIM + g J (l)2,l{x)cfx (1.1) 

where Atim represents the action of the c = 0.7 minimal CFT and g is a coupling constant, of 
dimension (mass)^, which describes the strength of the perturbation. It was shown in [|l| that 
the action ( |1.1|) describes a completely integrable QFT in the sense that it possesses an infinite 
number of non-trivial local integrals of motion. 

This model has been already studied by several authors Their results (with some vi- 

olations in chronological order) are summarised below. The truncated conformal space approach 
(TCSA) calculations by Lassig, Mussardo and Cardy showed that the theory (|1 . 1| ) has two 
degenerate vacuum states (labelled here by and 1) which correspond to the minima of the 
asymmetric double well Landau- Ginsburg potential. It was suggested that the particle spec- 
trum consists of a triplet of fundamental kink state^ l-^oi); l-^io)? aiid \Kqq) of the same mass 
m, which interpolate between these vacuum states (e.g. the kink \Kqi) interpolates between the 
vacuum states and 1, etc.). The fourth possible kink state \Kii) is not present due to the ab- 
sence of the Z2-symmetry between the vacuum states which is broken by the perturbation term 
in (|1 . 1| ) . This leads to a restriction on the allowable adjacent vacuum states, which (as noticed 
by Zamolodchikov in 0) is equivalent to the restriction imposed on adjacent spin states in the 
solvable lattice hard hexagon model (HHM) Therefore, assuming the above particle spec- 
trum is exact, one can construct the kink-kink ^-matrix from the associated Boltzmann weights 
of the (critical) HHM, since they satisfy the same Yang-Baxter equation. The normalisation 
of the S'-matrix is then determined by imposing unitarity, crossing symmetry and bootstrap 

requirements. The resulting S'-matrix is given explicitly in the next Section. Alternatively, (and, 
^ Note that the vacuum labels and 1 of paper || have been interchanged here for consistency with the 
usual convention of spin labelling in the HHM jT^ ]. 
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in fact, earlier) Smirnov |^ has conjectured S-matrices for all 0i^2 and 02,1 perturbed minimal 
models expressing them in terms of certain (appropriately normalised) RSOS projections of the 
^2^'' R-matrix of the Izergin-Korepin model . In particular, for the special case of the PTIM 
we are interested in here his result leads exactly to the same particle structure and S-matrix 
as discussed above. It must, however, be emphasised that both of these approaches are still 
based on conjectures and as such require further verification. 



In this paper we apply the thermodynamic Bethe Ansatz (TBA) approach [|TT],|T2| to the 
conjectured ^-matrix of the PTIM and calculate the ground state energy E{R) in the finite- 
size geometry with the spatial coordinate compactified on a circle of circumference R. The 
calculations utilise the above relationship between the PTIM and the lattice hard hexagon model 
allowing the use of known results from the analytic Bethe- Ansatz solution of the latter. The 
ground state energy E{R) is expressed in terms of the solution of the TBA integral equation of 
an apparently new type. We show that in the "ultraviolet" limit R —y the ground state energy 
scales as 

E(R) ~ (1.2) 

with c = 0.7, exactly as one expects from (|1 . 1|) . We also solve the TBA equations and evaluate 
the ground state energy E{R) numerically. The results are in a good agreement with those 
obtained by the TCSA [^,^ and from the Conformal Perturbation Theory (CPT) 



2 The Bethe- Yang equations 



As mentioned above the kink-kink S-matrix in the PTIM is expressed in terms of the Boltzmann 
weights of the critical HHM. The latter is an "interaction-round-a-face" model on the square 
lattice. With a suitable normalisation convenient for our purposes the Boltzmann weights of 
this model can be written in the form [1^ 
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where u is the spectral parameter and fi = tt/5. Note that this model is also solvable for other 
values of fi such that S/u/vr = integer, where it corresponds to some non-unitary QFT (see Sect. 
6 below). For the above normalisation of the Boltzmann weights the partition function per site 
of the model in the limit of the infinite lattice is given by |T3[ 

sinh(7r — fi)x sinh 2ux 



k{u, jj) = exp 



In particular for = 7r/5 it reads 



2x sinh ttx cosh fix 



dx 



< Reu < fi. 



(2.6) 



k{u, 7r/5) 



sin(7r/5 + u) sin(27r/5 — u) 



sin(7r/5 — u) sin(27r/5 + u) 
The kink-kink S-matrix proposed by Smirnov ||^ can then be written as0 



(2.7) 
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where 6 is the rapidity variable, po = 1; Pi = 2cos;U, while the other parameters are 



TT 
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(2.9) 



The normalisation factor in ( |2.8| ) is given by 

R{e) = ^{Xd,fi)-'FccD{0) 
where Fcdd{0) denotes the CDD factor 



(2.10) 



FcDD{0) = -F_.{d)F.{d), F^{9) 



sinh 6 + i sin an 
sinh 6 — i sin air 



(2.11) 



^ We follow p for the spin labelling convention in the elements of the S'-matrix, except for the aforementioned 
interchange of vacuum labels and 1. 



4 



It was shown @,|^,|^ that the ^'-matrix ( p.8| ) satisfies all the unitarity, crossing symmetry and 
bootstrap requirements. 

Consider the state of kinks distributed along a large spatial circle of the length L so that 
their average mutual distances are much greater than the correlation length ^ = 1/m of the 
system, where m is the kink mass. Since the kink scattering processes are of the factorised type, 



this state can be described in terms of the Bethe wave function |]rT|, [|T6| 



^(^1, ....On) =J2^^''^\Ka,aM)Ka,a,{02) ■ ■ • K^^aA^N)) (2.12) 

where ■ip^"'^ denotes the "colour" wave function depending on the set of interkink vacuum labels 
{a} = {«!, a2, . . . , «Ar}, tti = 0, 1, along the circle. Note that the absence of the kink Ku in the 
theory implies that the sum in (|2l2|) is taken over only those sequences of {a} which do not 
contain two consecutive a's equal to 1 in any place. 

Consistency of the imposed periodic boundary conditions requires that 

e'^^'^''^'''' E nOk, ^1, . . . , ^7v){S^^"^ = k = l,...,N (2.13) 

{a} 

where 

T{e; e,,..., ^^)|f| = n s;:^{e - e,) (2.14) 

1=1 

can be interpreted as the transfer matrix of a two-dimensional lattice model acting in the space 
of colour wave functions ip^"'^ The equations ( p. 131 ) result in constraints on the set of the kink 



rapidities 61,62, ... ,0^ which are known as the Bethe- Yang (BY) equations 

g^mLsinhe,^^^^. e,^...^6N) = -l, k = 1, . . . , N (2.15) 

where A(^; 6i, . . . , 6n) are the eigenvalues of the transfer matrix (p.l4|) . From the definition (|2.8|) 
it is obvious that ( p.l4|) is just the transfer matrix of the (inhomogeneous) critical HHM. The 
eigenvalues of the latter were found in |]IB[ by solving the transfer matrix functional equation. 



Using this result one obtains 



where 



JV 

A{6- 6,,..., 6n) = AHHMiXe- 6,,..., 6n) H ^(^ " ^3) (2-16) 
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N -TT siniu — X6j) ^, , ^ . / /i , \ 

^'">-n .in(.-A^,-,) - 0(")=n-("-|-A".) (2.1V) 



and where iV is assumed to be even, and the trivial difference in the normahsation of the transfer 
matrix ( |2.14| ) with that of ref. [T3\ has been taken into account. The (complex) numbers ak, 



k = 1, . . . , N/2, are determined by the following Bethe-Ansatz equations (BAE) 



N n 

uJ'^ n S^i^k - 0,) = - n S2iak - ai) (2.18) 



where 
and 



n = - (2.19) 



' sm{Xe-fij/2) ^ ' 

Note that the equations ( 2.18| ) can be considered as the BAE for the n = N/2 up-arrow sector of 



the inhomogeneous six-vertex model [|l^ with "twisted" boundary conditions. The factor uj~^ in 
the RHS is then interpreted as the "twisting factor" determined by the value of the (horizontal) 
field applied to the six-vertex model. 

3 The Thermodynamic Limit 



According to the standard approach of the Thermodynamic Bethe Ansatz (e.g. ||T2[, [T^) the BY 



equations ( p.l5| ) allow one to determine the spectral density of states of the infinite size system 



and hence to compute the free energy fo{mR) of this system at a finite temperature 1/R. This 
free energy is then reinterpreted as the Casimir part of the ground state energy E{R) of the 
finite-size system on a circle of the circumference R 

E{R) - (bulk term) = Rfo{mR) (3.1) 

where fo{mR) is the infinite-size free energy per unit length and the linear in R bulk energy 
term is known explicitly (see eq. ( ^4.5| ), below). 

In the thermodynamic limit the number of kinks and the system size L simultaneously 
approach infinity. From earlier studies of the six vertex model it is well known |jl9[ that for 
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N ^ oo the roots which solve the BAE ( p.lSD approach certain asymptotic patterns in the 
complex plane which can viewed as collections of "strings". Each string is a set of roots with 
the same real part, symmetric with respect to the real axis of a and equally spaced along the 
imaginary axis. The number of roots in a string is called the "length" of the string. In general 
the string spectrum of ( p. 18 ) may be very complicated |19| since it heavily depends on the 
arithmetic properties of the constant fi. For the case under consideration the constant fi is given 
by the simple fraction yU = tt/S and the string spectrum of ( p.l8[ ) consists of only five different 
strings. These are the strings of lengths i = 1,2,3,4 (to be referred as 1-string, 2-strings, etc.) 
which are centred at the real axis a = and the "shifted 1-string" which is located on the 
line Qma = 7r/(2A) (mod vr/A) and shifted from the real axis by half of the period of the 
function Q{\6). With these strings the roots for a typical solution of ( p.l8|) are given by 

«i!i = /^f + ^(^+l-2H + '^(^) (modvr/A), m = l,...,i, j = l,...,m (3.2) 
and 

4"^=/^r + ^ (modTT/A), j = l,...,n_ (3.3) 

where f3j^^ and are (real) centres of the strings while n^, i = 1,2,3,4 and n_ are the total 
numbers of strings of these types. The correction term S{N) in ( p.2|) describes deviations of the 
roots from exact string positions. For large it vanishes pO|3 as 0(logA^/A^). In the case of 



the HHM we only require those solutions of ( |2.18|) where the total number of roots is equal to 
n = N/2 and hence 

^ TV 

Y^in, + n.= (3.4) 

e=i ^ 

For a generic value of the twist factor uj on the unit circle, |a;| = 1, all the above strings are 
thermodynamically significant in the sense that for typical solutions of ( |2.18| ) with n = N/2 the 
ratios rii/N, £ = 1, ... ,4 and n_/N remain finite as ^ oo. However, when uj approaches 
the specific values logtu = ±2/c7r/5, k = 1,2,3,4 as it is required in ( p.l7| ) and ( p.l8| ) the 
situation drastically changes. It turns out in this case that, for many of these solutions, a finite 
number of strings moves away to infinity and the condition (|3.4| ) can no longer be satisfied. 
More precisely, on the basis of numerical calculations it was conjectured [O] that the solutions 



■^Actually, the behaviour of S{N) for large N is more compHcated The asymptotics S{N) ~ 0(\ogN/N) 
quoted above are valid for the region Re|a| <C logiV used for the thermodymic limit below. 
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of p.l8| ), where the fractions of the 4-strings Ui/N and/or the shifted 1-strings n_/N are not 
vanishing as iV — oo, do not match the condition (|3.4| ) and therefore have no relevance to the 
HHM. This leaves only three thermodynamically significant strings of lengths £=1,2, 3. Their 
total numbers are restricted by the relation 

Y^ine = - + o{N). (3.5) 

£=1 ^ 

where the o(A^) term grows slower than for oo. Assuming that the roots take their 

limiting string form (|3.2|) one can rewrite the equation ( 2.15| ) as 

N 3 rii 

^6^"^^^^°'^^'= n ^(^^ - ^^o n n seiOk - ^f) = -i + o(i), (3.6) 

i=i 1=1 3=1 



where is defined in (p. 201) and the correction term o(l) vanishes as A^ ^ oo. The last 
equation looks very similar to the BY equation for the diagonal scattering theory. In fact, it is 
very convenient to formally interpret the string centres j3f^ appearing in ( p.6|) as rapidities of 
some quasi-particles. These quasi-particles do not correspond to any observable particles in the 
asymptotic scattering states but rather reflect the nontrivial spin structure of the kink scattering. 
In the thermodynamic limit the rapidities of kinks and quasi-particles form dense distributions 
which can be described in terms of continuous densities. Let Pq{9) denote the rapidity density of 
kink states and pi{[3) denote such densities for the quasi-particles. We normalise these densities 
by the conditions 

Po{e)de=N- L pe{l3)d/3=ni, i = 1,2,3. (3.7) 

-oo J —oo 

Following the standard procedure of the TBA one can now rewrite ( p.6|) as an integral equation 

— coshe = Po(^) + Po(^)-E / '^ok{0-e')pk{e')de' (3.8) 



where, as usual, po{0) denotes the density of the "holes" in the kink rapidity distribution 
and 

*oo(^) = ^delogR{e), = ^,o{e) = ^de\ogS,{9), j = 1,2,3. (3.9) 

Similarly the BA equations ( |2.18|) lead to, 

3 nOO 

= p^{e) + p,{9) - E / ^,k{9 - 9')pk{9')d9' (3.10) 
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where Pj{0), j = 1, 2, 3, denotes the densities of holes in the quasiparticle rapidity distributions, 



\l/jo(^) is given above in (|3^ ) and 



^.fcW = -TT-de E ^og[S,^k+2U0)Sj+k-2mm, J,k = 1,2,3. (3.11) 
27r« m=i 



The equations ( p.8| ) and ( p.lO| ) can be notably simplified. In fact, integrating the equation (|3.10|) 



with j = 3 over 9 from — cxd to cxd and using ( |3.7| ) and the completeness relation ( |3.5D one obtains 



PsiO) = 0. (3.12) 

This means that the corresponding density of holes identically vanishes, p^iO) = 0. Therefore 
the equation ( ^.1(J| ) with j = 3 can be used to exclude the density p^i^O) from ( |3.8| ) and ( |3.10| ) 



expressing them through the remaining densities Pj{0), j = 0,1,2. Another transformation, 
especially useful for a simplification of the integral kernel in these equations, consists in the 
following. The densities Pj{0) and Pj{6) enter the integral equation ( ^.l(j| ) in a non-symmetric 
way. Indeed, the hole densities pj{0) enter only the free term in (|3.10| ) while the "particle" 



densities Pj{0) enter the free term and also the integral term in ( p. 101 ). This arrangement can 



be changed. In particular, multiplying the equations (|3.8|) and ( |3.10 ) with a suitable matrix 



integral operator one can bring them to an equivalent form where the "particle" densities Pj{0), 
for j = 1,2 appear only in the free terms (the corresponding hole densities will then enter the 
free and the integral terms). Performing all the transformations described above and redenoting 
the set of densities, 

ao{9) = po{9), MO) = PoiOy, a,{9) = p,{9), ^9) = p,{9), j = 1,2, 
one obtains 

6j,o cosh^ = aj{9) + aj{9) + sj ^ / '^jA^ - 0') (^k{9') d9', j = 0, 1, 2. (3.13) 

7, nJ~O0 



m 
2tt' 



k=0 

The new kernel reads 



= + Sj,k~i)MO) - So SjflSkfiMO) (3.14) 



^o(^) = ^ U V '^i^^) = 7^de\ogFcnD{0) (3.15) 

2p cosh(7r|A|6'//i) 2m 
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where Fcdd{0) is given by and the quantities Sj, j = 0, 1, 2 are sign factors 



.0 = s, = s, = 1. (3.16) 

Of course, in the case considered where A is given by (|2l9|), the value of Sq is fixed to be Sq = — 1. 
We have given the formula p.l6| ) just to illuminate the origin of these sign factors in p.l3| ). 



Thus, we have shown that in the thermodynamic limit the scattering state of the system of 
kinks in the PTIM is described by the three pairs of the rapidity and the hole densities (Tj{6), 
&j{0), j = 0, 1,2, which are constrained by the integral equation ( |3.13| ). In the next section we 
will use this result to compute the equilibrium free energy of the system at a finite temperature. 

4 The ground state energy of the finite size system 

In the thermodynamic limit, the free energy for the system of kinks considered in the previous 
section is given by the functional 

j^[a,{e),a,{e)] = s[a,{e)] - Ls[a,{e),a,{d)], (4.i) 

where 1/R acts as the temperature parameter 

/oo 
cosh9ao{9)de, (4.2) 
-oo 



is the energy of the system and S is the combinatorial entropy |^1| for the given set of the 
"particle" and "hole" densities aj{9),aj{9), j = 0, 1,2. 

Following the standard calculations the equilibrium free energy 

foimR) = - = -—— / cosh 9 log (1 + e~'''^^^ ) d9 (4.3) 

is obtained by minimising the free energy functional ([4.1D with respect to the above set of 



densities. The pseudo-energies 

e,{9)=\ogia,{9)/a,i9)) 
are determined by the TBA integral equations 



£■(9) = Sjfi mR cosh 9 + J2sk / <I>j,fc(^ - 9') log (1 + e-'"^^ ^ ) d9' (4.4) 



k=0 
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where the kernel ^jk{0) and the sign factors are given by ( p.l4[ ) and ( p.l6| ) respectively. 
Combining (|4.3| ) with (|3.1j ) and inserting the known bulk energy term, one obtains the ground 
state energy E{R) of the finite size system defined on a circle of circumference R 



E{R) =eR-— cosh ^ log (1 + e'^''^^) ) dO . 



(4.5) 



where [14 



24 sin(^: 



-0.41560346108. . .m^ 



(4.6) 



As usual, the leading asymptotics of E{R) in the ultraviolet limit, R 0, can be calculated 
by using the well known "dilogarithm trick" ||22|. In this limit the pseudo-energies ej{6) tend to 
finite constants ej{0) inside the region — In (^:^^ ^ 6* ^ In These constants satisfy the 

set of algebraic equations obtained by setting R = in ( [4.4|) . Bringing them into a symmetric 
form one gets 



y 



(0) 



k=0 



(o)yifc 



(4.7) 



where we have defined yl^^ = exp(— £0(0)) and y^'^' = exp{ej{0)), j = 1,2. The coefficients / 
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2 

3 

1 
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The equations (^77| ) have the unique positive solution 

V2, = 1, 



(0) 

yo 



1 



(4.9) 



Next, in the limit \9\ 00, the pseudo-energy eo(^) ~^ 00 while the other two pseudo-energies 
£1(6'), £2(6') tend to finite constants. These constants satisfy the set of equations 



(00) 

yl 



(00) 

y2 



1 + {yD 



(oo)x-l 



(4.10) 



where we have defined y'f^'' = exp{ej{oo)). These equations have the unique positive solution 



yir^ = 0, 



(oo) 

yi 



(oo) 

y2 



V5-1 



(4.11) 



where for convenience we have defined yQ°°^ = exp(— £0(00)) = 



11 



Following the standard calculations (see e.g. [|T^) one can deduce the leading i? — > asymp- 
totics of the ground state energy 

2 



3=0 



where 



1 

2 Jo 



C (fiyf : 



\n{t) ^ ln(l - t) 



dt 



(4.12) 



(4.13) 



1-t ' t 

is the Rogers dilogarithm function and f{x) = 1/(1 + x). With the simple dilogarithm identities 



C{x)+C{l-x) 
the equation ( [4.12| ) reduces to 



C(l] 



C(l] 



To 



E{R) 



77r 



60R 



(4.14) 



(4.15) 



exactly as one expects from ( |1.2| ) with the central charge c = 0.7 of the tri-critical Ising CFT. 
The following few terms of the short distance expansion of E{R) given by ( |4.5| ) can easily be 
determined numerically. This is done in the next Section where we compare the TBA results 
with the those of the Conformal Perturbation Theory (CPT)0. 



5 Numerical results 



The TBA equations ( [4.4|) can be solved numerically using a simple iterative procedure. The cor- 
responding ground state energy E{R) is then determined by numerically evaluating the integral 
in ( |4.3| ). Some resulting values of E{R) in the region 10^^ < mR < 10 are given in Table 1 
where they are compared to the corresponding values obtained from results of the Truncated 
Conformal Space Approach (TCSA) |2^. 

The TCSA calculations for the PTIM were performed in |^,^], however the required numer- 
ical data was omitted there. We are indebted to Prof. G.Mussardo who provided us with the 
computer program [2^] which generates the TCSA data presented in Table 1. Note that to relate 
the TCSA data back to TBA results one has to use the exact relation |T^ between the kink 



It should be noted that bulk term in (4.5) exactly compensates the next-to-leading, linear in R, term in the 



short distance expansion of the integral in (4.5). In fact, the main correction term to (4.15) is of order 0(i?^/^), 
exactly as predicted by CPT. 
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mR 


E{R)/m (TBA) 


(TCSA) 


(CPT) 


0.00001 


-36651.91429205 


-36651.914292 


-36651.9142920555 


0.000025 


-14660.7657173015 


-14660.765717 


-14660.7657173015 


0.00005 


-7330.38285968217 


-7330.382860 


-7330.38285968217 


0.000075 


-4886.92190775210 


-4886.921907 


-4886.92190775210 


0.0001 


-3665.19143229428 


-3665.191432 


-3665.19143229428 


0.00025 


-1466.07658143979 


-1466.076578 


-1466.07658143979 


0.0005 


-733.038309061764 


-733.038300 


-733.038309061764 


0.00075 


-488.692229111026 


-488.692214 


-488.692229111026 


0.001 


-366.519198155448 


-366.519177 


-366.519198155448 


0.0025 


-146.607830808513 


-146.607765 


-146.607830808515 


0.005 


-73.3042415739432 


-73.304084 


-73.3042415739653 


0.0075 


-48.8699046289368 


-48.869643 


-48.8699046290281 


0.01 


-36.6528965534080 


-36.652522 


-36.6528965536577 


0.025 


-14.6638535325716 


-14.662677 


-14.6638535387397 


0.05 


-7.33772690834948 


-7.334929 


-7.33772697813054 


0.075 


-4.89911302398585 


-4.894469 


-4.89911331240653 


0.1 


-3.68265799877659 


-3.676005 


-3.68265878811873 


0.25 


-1.52096720347295 


-1.500051 


-1.52098667429440 


0.5 


-0.86341841907386 


-0.813667 


-0.86363725527810 


0.75 


-0.70459594978559 


-0.621994 


-0.70548946787533 


1.0 


-0.67473619457720 


-0.556359 


-0.67713759026010 


2.5 


-1.07537196387997 


-0.702079 


-1.12306270607538 


5.0 


-2.08009259839356 


-1.182952 


-2.39571841186687 


10.0 


-4.15604421536308 


-1.960064 


-5.56031580714459 



Table 1: Some numerical values for the ground state energy E{R) (given in units of mass m) 
computed from the TBA equations ( [4.4| ),( |^75| ), from the TCSA method 1^^, and from the first 
three terms of the CPT short distance expansion ( |5.2| ), for different values of the dimensionless 
variable mR. 



mass m and the coupling constant g in (|1 . 1|) 

2W4/3r(5/9) / 2r(ll/16) 



m 



^ 4.9277912244. . .c/ 



8/9 



(5.1) 



r(2/3)r(8/9) vr(5/i6)r(3/4) ; 

The difference between the TCSA and TBA results given in Table 1 (seen more clearly in 
Fig. 1) is attributed to both the inherent inaccuracy of the TCSA truncation procedure (used 
to approximate the Hilbert space of energy states) which increases with R, and also to the fact 
that the dimension A = ^ is very close to the critical value of A = 1, above which divergence 
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-1.5 



DC 
LU 




-2.5 



Figure 1: The ground state energy E{R) (given in units of mass m) computed from the TBA 
equations ( [4 .41 ), ([4 .51) , from the TCSA method |2^, and from the first three terms of the CPT 
short distance expansion ( |5.2| ), plotted against the dimensionless variable mR. 



problems are known to occur. It is the latter problem which is thought to be responsible for the 
incorrect value of the bulk linear slope in the TCSA data at large R (as exemplified in Figure 1). 
We refer the reader to the original paper where these problems are more thoroughly discussed. 

In Table 1 (and Fig. 1) we also compare our numerical results for E{R) with the first three 
terms of the short distance expansion 



TT 



- -6R 



10 V27r 



(5.2) 



obtained from Conformal Perturbation Theory (CPT) [12, 16, 25|. In our case the coefficients 62 
and 63 read explicitly |14 



12 



(^) r (i) 



135.9255988883... 



0. 



(5.3) 
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Considering these two coefficients as fixed and using the TBA results it is not difficult to make 
a numerical ffi for a few following coefficients in (^.2|) . In particular, we found that 

64^-250.3..., 65- -0.58..., 66 ^2.3... X lOl (5.4) 



6 Discussion 

6.1 There are several other models of integrable QFT where the particle spectrum only consists 
of kink excitations with S'-matrices related to the Boltzmann weights of the critical HHM. The 
most known among them is the RS0S(4) scattering theory arising from the perturbation of 
the (c=0.8) minimal CFT A^s/e by the operator ^13. The S'-matrix of the general RS0S(j9) 
scattering theory, p = 3, 4, . . . , just coincide with the "unitarised" Boltzmann weights^ of the 
critical lattice RSOS model of Andrews, Baxter and Forrester (ABF) I^G] . This is an interaction- 
round-a-face model on the square lattice where admissible neighbouring spin states belong to 
adjacent sites of the Ap incidence diagram. The transfer matrix of this model commutes with the 
"spin reversal" operator generated by the reflection symmetry transformation of the incidence 
diagram and therefore it splits into a direct sum of two transfer matrices acting diagonally in 
the even and odd (with respect to the above symmetry) subspaces of the space of states. In the 



case p = A each of these two transfer matrices coincide with the transfer matrix of the hard 
hexagon model (modulo minor differences in boundary conditions which are irrelevant in the 
thermodynamic limit). Therefore, for the purposes of the TBA calculations, the kink S'-matrix 
of the RSOS (4) model can be replaced by an appropriate "hard hexagon" type S'-matrix which 
IS given by the formulae (l23),(iD and (|2lO| ) with A and FcDD{d) replaced by 



/i = |, ^ = -^' FcDD{e) = l, (6.1) 

and all other quantities to remain unchanged. The calculation of Sect. 2.4 requires no modifica- 
tion for this case, one just needs to specialise the final results ( [4.4|) and ( |4.5|) . 
The ground state energy is 

E{R) = (bulk term) / cosh 9 log (1 + e-""^^) ) d9 . (6.2) 

2ti J-oo 



^The Boltzmann weights should also be supplied with simple gauge transformation factors which ensure the 
crossing-symmetry of the 5-matrix. 
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where the pseudo-energy eo{9) is determined by the following TBA integral equations 
ej{9) = SjflmRcosh9 + Y, ^j,k{0 - 9') \og{l + e-"^^^ ^ ) d9' , j = 0,1, . . . ,p - 2. (6.3) 

1. n J —C>0 



with p = A. The kernel $j,fc(^) is given by the same formula ( |3.14|) which, with account of ( |6.1| ), 
reads exphcitly 



Note that the sign factors sj which are present in (|4.4| ) do not appear in (|6.3|) because the A in 
( |6.1[ ) has the opposite sign to the one in (|2.9|). 

The main technical point of our calculations was the diagonalization of the transfer matrix 
( p.l4| ) by using the analytic Bethe-Ansatz solution [|13| for the HHM. In fact, the aforementioned 
paper contains the results of this diagonalisation procedure for the complete hierarchy of ABF 



RSOS lattice models |2g] (and their higher fused generalisations |23|). The calculations of the 



present paper are also easily extended to the case when the kink S-matrices are expressed 



through the Boltzmann weights of these more general models |28|. In particular, for the general 



RSOS(p) scattering theories one gets in this way the TBA equations ( p73D , which are identical 



to those conjectured by Zamolodchikov ||I6[. Note also, that essentially the same set of TBA 

equations (or more precisely their lattice counterparts, which differ from ( |6.3| ) only in the form 

of the energy term mi? cosh 6') arose previously [l^ in the analysis of finite size corrections to 

the critical ABF RSOS models. 

6.2 Another integrable QFT where the scattering theory is related to the HHM is obtained 

from the c = —3/5 minimal CFT Ais/^ perturbed by the operator 02, 1 with dimension A = 3/4. 

This QFT is also covered by Smirnov's conjecture ||^. It was specialised to this case in p9[] , 

resulting in a kink S'-matrix that can be written exactly as in (2.10) and (2.12) with 

37r 3i sinh^ + 2sin(7r/3) 

= ^, A = - — , Fcdd{0) = — r-7 . . , (6.5) 

5 5 smliy — zsm(7r/3j 

Evaluating the integral in ( p.6|) for /i = 37r/5 and substituting the result in the normalisation 



factor (p.lO|) one gets explicitly 



R{9) = >^ '-4 \^ (6.6) 

^ ^ sin + f ) sin (f + f ) 

The calculations of Sect. 3 and Sect. 4 can be repeated in this case with certain modifications. 
The main difference is related to the fact that for fi = 37r/5 the string structure of the solutions 

16 



of the BAE equations ( |2.18| ) changes. Numerical calculations suggest that in this case only 
1-strings in (3.2) and the shifted 1-strings (3.3) are thermodynamically significant. We believe 
this statement is correct and claim it as a conjecture. The rest of the calculations closely parallel 
those given in the case /i = vr/S and requires no further assumptions. As a result one obtains 
that the ground state energy E{R) is still given by the expression (3^) while the pseudo-energy 
eo{9) is determined by the following TBA equations 

So{e) = mR cosh 6+ H <^{e - 6') (log (1 + e-'^^'^'^ ) - log (1 + e'^^^^') )) dO', 

ei{e) = + r - 6') (log (1 + e~'^^''^ ) - log (1 + e-^°(^') )) dO', (6.7) 

J —oo ^ ' 



where the kernel is given by 



MO) = (6i 



TT sinh(36') 
Note also that 

^{e) = ^de\ogFcDD{e) (6.9) 

ZTT 

where Fcdd is the CDD factor defined in (|6.5D. 

Apparently the same TBA equations for this QFT were conjectured in from different 
arguments. 

6.3 Yet another related QFT is obtained as the perturbation of the minimal CFT A^s/e by the 
operator 01^2 of dimension A = |. Again the 5- matrix of the fundamental kink states (as given 
by Smirnov's conjecture Q]) is expressed through the Boltzmann weights of the HHM (2.3-2.7) 
with fi = 7t/5. However, unlike all of the above examples, the bound state state structure of this 
QFT appears to be extremely complicated. Apart from the fundamental kinks it contains higher 
kinks and breathers as bound states of the fundamental kinks. These higher kinks and breathers 
produce even more new kink and breather states and so on. A few of the lower kink and breather 



states were found in |Q . It would be interesting to see whether this bootstrap procedure closes 



^In Ig^l the TBA equations for the -M^/^ minimal CFT perturbed by the operator 02, i were obtained by 
"orbifolding" of the TBA equations associated with the product of A2 Cg) A2 Dynkin diagrams. The latter were 
previously conjectured in |3^,^ (in a rather general form) as field theory counterparts for the TBA equations of 
lattice RSOS models connected with the simply laced algebras obtained in ref. [Q. It should be noted that the 
relevant equations in refs. ]3^-|3^ contain some consistent misprints. However, tracing their derivations back to 



the original reference |3^] one can recover the TBA equations (3.7) exactly as given here. We thank F. Ravanini 
and R. Tateo for explaining this issue. 
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at all. Particularly, some preliminary estimates show that the spectrum of this theory contains 
at least a dozen particles. In any case all of the kink-kink S'-matrices are expressed in terms of 
the HHM Boltzmann weights (2.3-2.7) and therefore the TBA calculation of the present paper 
could be generalised to this case provided the bootstrap program is completed. 



The above considerations were restricted only to the calculation of the ground state energy 
in the finite-size IQFT's by the TBA method. Recently, a few new approaches were introduced 



35 -39 1 which allow one to calculate also excited state energies. In ref. |37| it was shown that 



the functional equations for the IQFT "commuting transfer-matrices" introduced in |^ can 
be transformed [0 to the integral equations which generalise the equations of TBA to the 
excited states. The results of apply to massive IQFT's obtained by 0i_3 perturbations of the 
minimal CFT's which are related to the Uq{sl{2)) quantum algebra. A suitable generalisation 
of this approach to treat the 0i_2 or 02,1 perturbations of the minimal CFT's (in particular, the 
PTIM considered above) should be associated with the g-deformed twisted Kac-Moody algebra 
Some results in this direction were obtained in We hope to address this problem in 
future publications. 
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